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It has been argued that the small perturbations to the homogeneous and isotropic configurations 
of a canonical scalar field in an expanding universe do not grow. We show that this is not true 
in general, and clarify the root of the misunderstanding. We revisit a simple model in which the 
zero-mode of a free scalar field oscillates with high frequency around the minimum of the potential. 
Under this assumption the linear perturbations grow like those in the standard cold dark matter 
scenario, but with a Jeans length at the scale of the Compton wavelength of the scalar particle. 
Contrary to previous analyses in the literature our results do not rely on time-averages and/or fluid 
identifications, and instead we solve both analytically (in terms of a well-defined series expansion) 
and numerically the linearized Einstein-Klein-Gordon system. Also, we use gauge-invariant fields, 
which makes the physical analysis more transparent and simplifies the comparison with previous 
works carried out in different gauges. As a byproduct of this study we identify a time-dependent 
modulation of the different physical quantities associated to the background as well as the pertur¬ 
bations with potential observational consequences in dark matter models. 

PACS numbers: 98.80.-k, 98.80.Jk, 95.35.+d, 03.50.-z 


I. INTRODUCTION 

It is difficult to overestimate the relevance of the Jeans 
instability in modern physical cosmology. In order to un¬ 
derstand the emergence of cosmic structure, we need a 
mechanism that transforms the nearly homogeneous and 
isotropic early universe we infer from e.g. the cosmic mi¬ 
crowave background observations, to the highly clumped 
one we can see today at the scale of galaxy clusters and 
below. In the standard cosmological scenario this transi¬ 
tion is possible thanks to the instability of a scalar mode 
that appears when we couple dark matter (DM) to grav¬ 
ity, the gravitational Jeans instability. 

In a universe dominated by a barotropic perfect fluid 
with equation of state p = p(e) -C e, the behavior of 
the small perturbations in the energy density depends 
crucially on the speed of sound, c 2 = dp/de, and the 
Hubble radius, H~ 1 (throughout this paper we use nat¬ 
ural units such that c = h = 1). Roughly speak¬ 
ing, we can distinguish three different regions in Fourier 
space [lj: i) On scales smaller than the Jeans length, 
Xj = c s (7r/Ge 0 ) 1/2 ~ c s H ~ l , the contrast in the energy 
density oscillates with damping amplitude, due to the 
stabilizing effect of pressure and the expansion of the uni¬ 
verse, respectively. Here G is the Newton constant and £q 
is the background matter density, ii) Above this length 
scale but still below the Hubble radius, self-gravity dom¬ 
inates and the contrast in the energy density grows: the 
Jeans instability comes into play. Hi) Finally, at scales 
larger than the Hubble radius a relativistic understand¬ 


ing of the problem shows that the contrast in the energy 
density freezes. 1 

Perfect fluid structure formation then demands matter 
with a low speed of sound in order to have a window, 
c s /H < A < 1/H, where the Jeans mode is released and 
the perturbations can grow. Cold dark matter (CDM) 
represents the simplest realization of this scenario, for 
which one assumes c 2 ~ 0 for all relevant modes. 2 

For a canonical scalar field we have c 2 = 1, see e.g. 
Ref. Q, and it is usually argued that this implies that 
small perturbations in the energy density do not grow [|[. 
This would seem to be a very serious argument against 
the whole scalar held DM program Q, so we find it 
mandatory to clarify the issue. 

In this paper we will consider the simplest situation: 
a universe dominated by a real massive scalar held ip 
satisfying the Klein-Gordon equation, (□ + m 2 )ip = 0. 


1 Actually all these statements are gauge dependent. We can al¬ 
ways choose to work in e.g. the uniform density gauge, where the 
contrast in the energy density vanishes identically at all scales, or 
in e.g. the synchronous gauge, where the contrast in the energy 
density grows even for those modes larger than the Hubble ra¬ 
dius. It is only in terms of the conformal-Newtonian gauge that 
the behavior outlined in the previous paragraph makes sense, and 
it is only in this gauge that we can easily compare our results 
with those sketched in the three points above. 

2 If DM consists of collisionless particles it is more appropriate to 
talk about a free-streaming, rather than a Jeans, length; however, 
the idea is similar, see e.g. Sections 10.2 and 10.3 in Ref. Q] for 
a discussion. 
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FIG. 1: Behavior of the linear perturbations in the contrast 
to the energy density for different wavelengths A. The Jeans 
length A j depends on the matter content in the universe. For 
a perfect fluid A j ~ c s H -1 , and structure formation demands 
c 2 -C 1. However, when a canonical scalar field oscillates with 
high frequency around the minimum of the potential we have 
A j ~ m -1 , even though (? s = 1. As usual this picture should 
be understood in the conformal-Newtonian gauge. 

Here the box is the d’Alembert operator in four space- 
time dimensions, and m is the mass of the scalar par¬ 
ticle [this corresponds to a scalar field potential of the 
form V(ip) = m 2 <p 2 / 2, so that the mass is defined as 
to 2 = d 2 V/dip 2 \. The more interesting case with a com¬ 
plex field including the self-interactions and the presence 
of other matter components (e.g. radiation, baryons, etc) 
will be presented elsewhere. We will show that when the 
scalar field is slowly rolling down the potential, linear 
perturbations cannot grow, in accordance with common 
wisdom. However, when the scalar field oscillates with 
high frequency m H around the minimum of the po¬ 
tential term, the Jeans length is not given by the naive 
value c s H ~ 1 one would guess from a perfect fluid anal¬ 
ogy. Instead it is determined by the Compton wavelength 
of the scalar particle, c s m -1 , and the evolution of per¬ 
turbations larger than this scale almost mimics that of 
the standard CDM scenario even though c 2 = 1 (see Fig- 
ure|T|for details). The main difference with respect to the 
standard CDM evolution (apart from the appearance of 
a nonvanishing Jeans length), is the presence of a time- 
dependent modulation of the different physical quanti¬ 
ties. We will present below both an analytical argument 
based on expansions of the solution of the relevant cos¬ 
mological equations, and results from simple numerical 
simulations. 

Similar results have beenpresented before in the litera¬ 
ture, see for instance Refs. 0 (as far as we know the grav¬ 
itational instability of a canonical scalar field in the 
context of a static universe— was reported for the first 
time by Khlopov, Malomed and Zel’dovich in Ref. 0 )- 
However, to our knowledge, we present for the first time 
a description of the problem in terms of gauge-invariant 
fields. Furthermore, note that our analysis does not rely 
on time averages and/or fluid identifications. Instead 
we present a well-defined series expansion that makes it 
possible to find analytic solutions order by order in the 
expansion parameter: it is well known that the time av¬ 
erage of the product of two functions is not in general the 
product of the time averages of those functions, and also 
that a scalar field is not a perfect fluid 3, so one needs 
to take some care with the standard procedure. With 
this formalism the results emerge more naturally than in 
previous works, and it is convenient in order to identify 
the root of the misunderstanding; see Eq. m and the 
paragraphs below. 


We will follow the notation in Chapters 7 and 8 of 
Ref. [ij; in particular we use the signature (+,—,—,—) 
for the spacetime metric. We highly recommend this ref¬ 
erence to the reader interested in the details about some 
of the definitions and conventions below. One should also 
mention that recently the problem of structure formation 
with a massive scalar field has been studied by perform¬ 
ing full nonlinear numerical simulations of the Einstein- 
Klein-Gordon system, both in the relativistic Q and the 
nonrelativistic [10] regimes. However, we believe that by 
studying the problem from the point of view of a mode 
analysis in perturbation theory one can more clearly sep¬ 
arate the relevant physical mechanisms that come into 
play at different scales. 

II. THE HOMOGENEOUS AND ISOTROPIC 
BACKGROUND 

At very large scales the universe is (nearly) homoge¬ 
neous and isotropic; that makes it possible to introduce 
the idea of a homogeneous and isotropic background. 
According to the current cosmological observations this 
background can be described in terms of a flat Robertson- 
Walker (RW) metric of the form 

ds 2 = a 2 (dy 2 — dx 2 — dy 2 — dz 2 ) . (1) 

Here 77 is the conformal cosmological time, ( x , y , z) a 
spatial coordinate system comoving with the expansion, 
and a{rj) the scale factor. Conformal time 77 is related 
to the standard comoving cosmological time t through 
t = f adrj. The expansion rate is codified in the Hubble 
parameter, H = T~L/a = a'/a 2 , with the prime denot¬ 
ing the derivative with respect to conformal time. The 
spacetime symmetries in Eq. © guaranty that the back¬ 
ground field cannot depend on the spatial coordinates, 
so that ip(r],x) = (po(v )■ Under these assumptions the 
Klein-Gordon equation simplifies to 

¥>o + 271(^0 + a 2 m 2 ipo = 0. ( 2 ) 

There are two different timescales in Eq. ©: on the 
one hand that associated to the cosmological expansion, 
T~L -1 = (ai?) -1 , and on the other that defined by the 
mass of the scalar field, (am) -1 . In this paper we will 
concentrate on the case when the expansion rate of the 
universe is very slow when compared to the characteristic 
time of oscillation of the scalar field, that is (ai7) -1 
(am) -1 , which implies H <C m. This inequality is always 
satisfied at late times. We will briefly discuss the regime 
H m later on in this section. 

In order to proceed we propose a solution of the form 

Vo(v) = l Po(v)[ l Po l (v) + 0 2 (H/m)] , (3a) 

Utn) = U H {r,) [l + 0(H/m)\ , (3b) 

a(rj) = a H (77) [l + 0 2 (H/m)] . (3c) 

Here functions with a superscript H vary on cos¬ 
mological timescales, ( f H )' ~ aHf H , whereas those 
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with a superscript m vary on timescales given by the 
mass of the scalar field, ( f m )' ~ amf m . In princi¬ 
ple, higher-order terms oscillate in time with high fre¬ 
quency. Consequently, the derivative of one of those 
terms is not necessarily suppressed in the series expan¬ 
sion, e.g. d[0(H/m)\/dr] ~ aH. That is the reason for 
which the series expansion of the Hubble parameter H 
has a linear term in H/m , even though such a term is 
not present in the expression for the scale factor a. 

Introducing the ansatz © into the Klein-Gordon equa¬ 
tion we obtain 

<Po(v) = [sin(mi) + 0 2 (H/m )] , (4) 


with A an integration constant, Mpi = l/V&i tG the re¬ 
duced Planck mass, and where we have fixed an arbitrary 
phase to zero. For later convenience, and with no loss of 
generality, from now on we will also choose A 2 = 24. 
From the above expression we see that the scalar field 
oscillates around p o = 0 with constant frequency to in 
comoving time t (so that the frequency will increase in 
conformal time r) as the universe expands), and an am¬ 
plitude that decays as 1 /rf. Introducing this expression 
for po(rj) into the Friedmann equations, 

U~ = — a 2 e 0 , (5a) 

H' - T-L 2 = —47rGa 2 (eo +Po) , (5b) 


one finds 


n( v ) 

a(v) 


— 1-(—^ sin(2mt) + 0 2 (H/m) 

77 |_ 4 \m) 

{mrjj 2 [l + 0 2 (H/m )] . 


( 6 ) 

(7) 


For completeness we can also integrate the expression for 
the comoving time to obtain 


mt = / ma drj = 


(miff 


[1 + 0\H/m)}. ( 8 ) 


In the above equations we have used the fact that the 
background energy density and pressure are given by 


£0 


Po 


2 2 
-TO <p Q 


1 


P 0 = 


Po 


-m p 0 ■ (9) 


In particular we find that, to lowest order in the 
series expansion, the background energy density red- 
shifts with the inverse of the comoving volume, 
£0 ~ 1 /a 3 , whereas for the background pressure we ob¬ 
tain po ~ (1/a 3 ) cos(2 mt). Note that to this same order 
we cannot distinguish the expansion rate from that in a 
CDM universe, H = 2/ 77 , even tough |po| ~ £0 during 
the evolution. Interestingly, this does not depend on the 
mass of the scalar particle, as long as it is large enough 
when compared to the expansion rate of the universe. 

Since H/m ~ 2/( 777 . 77 ) 3 , the condition H -C to 
demands mrj 1. In terms of the scalar field, 



FIG. 2: Evolution of the conformal Hubble factor "H/m as 
a function of the conformal time mrj for a free scalar field 
oscillating around the minimum of the potential. We show 
for comparison the solution in standard CDM (solid line), 
and both the analytical approximation in Eq. JSJ (dashed 
line), and the full numerical solution of the Friedmann-Klein- 
Gordon system, Eqs. © and ©, (dotted line). In the figure 
the integration begins at rmj = lO 1 ^ 3 , with a = and 

H/m = 2 x 10~ 1/3 . 


Eq. (H above, we obtain |<^o| Mpi. That guaran¬ 
tees large values kjJ = |t?sr| = 2 Mp|/ip§ 1 of the slow- 
roll parameters [llj, with e sr = (Mpj/ 2){d v V/V) 2 and 
r) SI = M/ l (d 2 V/V), so we can safely conclude that the 
universe is not in a period of slow-roll inflation, as was 
evident from Eq. ©. 

At this point one should mention that the dominant 
terms in the solution of the Friedmann-Klein-Gordon sys¬ 
tem given by Eqs. Q, © and © coincide with those of 
the exact solution found in Ref. [l^l for a scalar field 
evolving in a universe dominated by a barotropic fluid 
such that H = 2/(3~/t), with 7 constant (when 7=1 
this background barotropic fluid can be associated with 
the average energy density and pressure of the scalar field 
itself). Something similar happens with the solution re¬ 
ported in Ref. (l3f | . obtained in terms of time averages. 
The main difference with our results is the fact that here 
we have made explicit the existence of higher-order terms 
in H/m , and in particular we have shown the first oscil¬ 
lating subdominant contribution to the Hubble parame¬ 
ter. Note however that this is not only a purely academic 
question: the inclusion of the higher-order terms will be 
crucial next for the understanding of the evolution of the 
small perturbations. 

Figure [2] shows the evolution of the conformal Hubble 
factor H/m as a function of the conformal time mrj for 
a free scalar field oscillating around the minimum of the 
potential. The solid line corresponds to standard CDM, 
while the dashed and dotted lines show the case of the 
scalar field using both the analytical approximation in 
Eq. ©, and a full numerical solution of the Klein-Gordon 
and Friedmann system, Eqs. m and ©, respectively [for 
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the numerical case we solve for T~L from Eq. and use 
Eq. (153 to monitor the numerical error]. Notice that the 
evolution for the case of the scalar field follows closely 
that of CDM with small oscillations around it. Initially 
we have H/m = 0.2 and the oscillations are still evident, 
but their amplitude decreases as time goes on. 

Next we will consider the behavior of the small pertur¬ 
bations around the homogeneous and isotropic solution 
in Eqs. 0, ©, and 0. In particular, we find that an 
unstable Jeans mode grows at scales between the Comp¬ 
ton length of the scalar particle and the Hubble radius 
of the universe; see Eq. CD below for details. 


III. THE LINEAR PERTURBATIONS 

In this paper we will not consider vector and tensor 
modes, and will concentrate on the scalar sector of the 
perturbations. After all this is the sector that contains 
the Jeans mode of the theory. With this assumption, 
the most general expression for the spacetime metric of 
a universe close to a flat homogeneous and isotropic RW 
one takes the form 


two different series expansions: one in the small pertur¬ 
bations of the spacetime metric and the scalar field, and 
the other in the expansion rate of the universe. We will 
soon introduce a new one in terms of the Jeans length.) 

In order to move forward we find it convenient to define 
the new quantities 

» = „(*?+$*), z = ^. (13) 


The gauge-invariant field v(ri,x ) is usually known as the 
Mukhanov-Sasaki variable, and the function z depends 
only on the spacetime background, z = z{rf). Note that, 
leaving the scale factor aside, the Mukhanov-Sasaki vari¬ 
able represents the scalar field perturbation evaluated 
in the spatially flat gauge ip = E = 0 [l4|. Whenever 
the scalar field dominates the evolution of the universe, 
Eqs. (1 12 all and (I12bl) simplify to 


*(?H^ 2 ©'- <“*> 

(ir)'-* 0 *-®- < 14b > 


ds 2 = a 2 {(1 + 2<p)dr/ 2 + 2dtB dx l dri 

— [(1 — 2 ip)5ij — 2 didjE\ dx l dxi} . (10) 

For the perturbed scalar field we will also write 
p = Po + Sp. Here <p, B 1 ip, E and dp are functions of 
the spacetime coordinates r/ and x, with dp -C po and 
<p, B, ip, E -C 1. 

All of these fields depend on the choice of coordi¬ 
nates used to write the line-element in Eq. m, and 
they could in principle be describing fictitious inhomo¬ 
geneities. For this reason it is sometimes convenient to 
work with gauge-invariant fields, such as pj 


^ - (1/a) la(B-E')]', 

(lla) 

ip + n(B-E'), 

(lib) 

dp-p' 0 (B-E'). 

(11c) 


Note that the above fields coincide with the amplitude 
of the metric and the scalar field perturbations in the 
conformal-Newtonian (also known in the literature as the 
longitudinal) coordinate system, for which B = E = 0. 

To linear order in the new field variables, the 00 and Of 
Einstein field equations, essentially the Hamiltonian and 
momentum constraints, take the following form: 

AT - iU (T' + H$) = 47rGa 2 & 

= 47tG \p' 0 dp' + a 2 m 2 podp — p'o&\ , (12a) 

T' + "HT = 4wGp' 0 dp , (12b) 

where A is the Laplace operator in flat space, and where 
Ss = d£—£' 0 (B—E') is the gauge-invariant density pertur¬ 
bation. For a scalar field, and to this order in the series 
expansion, there are no anisotropic stresses and the i ^ j 
field equations fix 4> = T. (Note that so far we are using 


These two equations can be combined to obtain 

z" 

v" — c 2 Av — — v = 0 , (15) 

with c 2 = 1 analogous to a “speed of sound” for a canon¬ 
ical scalar field [3j. The main reason for this identifica¬ 
tion is the similarity of Eq. m with that obtained for 
a barotropic perfect fluid with sound speed c s ; see e.g. 
Eq. (7.65) in Ref. (TJ. Notice that even though here we 
are only interested in the case of a scalar field with no 
self-interactions, so far the analysis of the perturbations 
is general and valid for an arbitrary potential. 

Equation (fl5l) above fixes a characteristic length scale 
given by c s a\z" I z\~ x l 2 (in this paper we will work with 
comoving wavenumbers k, but we will talk about phys¬ 
ical wavelengths A, with A = 2ira/k). This length pro¬ 
vides us with a ruler to discriminate between large and 
small scales in the scalar field perturbations. The correct 
estimation of this length scale for the different physical 
situations —and in particular the relative size it takes 
when compared to the Hubble radius— lies precisely at 
the heart of the usual misunderstanding. 

As we will find soon, see Eqs. (l20l) and (1211) below, the 
scale c s a\z" / z\~ x ! 2 determines the Jeans length of the 
scalar field. However, this quantity is not always related 
to the naive guess c s H~ x one would expect from a per¬ 
fect fluid analogy. In order to clarify this point, let us 
consider two representative cases. When the scalar field 
is rolling down to the minimum of the potential we have 
p' 0 ~ aHpo, which implies \z" / z\ ~ ( aH) 2 . This fixes the 
characteristic Jeans length to the Hubble radius. Short- 
wavelength perturbations (when compared to this length 
scale) oscillate in space and time as v ~ sin(c s fc77 + k ■ x), 
with c 2 = 1. The Jeans mode is then stabilized and the 
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perfect fluid analogy seems possible, i.e. small pertur¬ 
bations to the homogeneous and isotropic solutions do 
not grow. Consider, for instance, the case of the infla- 
ton during a slow-roll regime, or a quintessence field in 
the present universe, as particular realizations of this sce¬ 
nario. 

However, when the zero-mode of the scalar field is oscil¬ 
lating with high frequency m^$> H around the minimum 
of the potential, we have instead cp' 0 ~ amtpo, which re¬ 
sults in \z"/z\ ~ (am) 2 . The new length scale is well 
inside the Hubble radius, making possible the growth of 
perturbations with aH < k < cj 1 am (remember that for 
a canonical scalar field c 2 = 1). Now the naive estima¬ 
tion c s H~ 1 for the Jeans length has nothing to do with 
the correct one, c s m~ l . 

Let us consider in more detail this second scenario. 
Notice that, using the solutions in Eqs. 0, ©, and 0, 
the “mass” term in Eq. (Tl5l) takes the form 


= —a m 


1 + 6 ( — ) sin(2mt) + 0 2 (H/m) 
V m 


■ (16) 


As mentioned before, this is a purely background quan¬ 
tity. For practical reasons we will work with periodic 
boundary conditions over a box of comoving size L. We 
can always take the limit L — > oo at the end of the calcu¬ 
lations. Under this assumption, the general solution to 


J 


Eq. (1151) can be written in the form 

v (fo J2 [ C k v k(v)e ik ' S + c.c. 

MO 


(17) 


Here Cr are some dimensionless integration constants 
that label the different possible solutions, c.c. denotes 
complex conjugate, and the functions v^(rj) satisfy the 
equation 

v l + w \( r i ) v %= 0 > with = fc2 - — • ( 18 ) 


Note that the boundary conditions demand ki = ‘hrriil X, 
with rii = ±1,±2,... and i = 1,2,3. The zero-mode 
k = 0 is already included in the description of the space- 
time background, and for that reason it does not appear 
in Eq. (fT71) . 

According to the Eq. (1181) the behavior of the mode- 
functions vArj) depends crucially on the relative value 
between the square of the wavenumber k 2 , and the back¬ 
ground function z" / z. For those modes smaller than the 
Compton wavelength, k am, we obtain (remember 
that in the fast expansion regime the Compton wave¬ 
length is always well inside the Hubble radius, and then 
am aH) 


vt(rj) = 7 = (1 + 0 2 (H/m , am/k )) exp [— i (l + 0 2 (H/m, am/k )) krj\ , 

V2fc 


(19a) 


whereas for modes larger than this quantity, k -C am, we find 


M = 


Mk) f 

\/2m 1 


B(k)z - \B(k)a ( — ) ( -77 
am J \ aH 


3 [ H 

sin(mt) — - ( — ) cos 3 (mt) + 0 2 (H/m) 


—iB x (k)z 


d(mrj) 




k 


2 \m 

3 ( H 

cos (mt) -( — ) sin 3 (mt) + 0 2 (H/m) 

2 V m 


+0 2 (H/m,k/am) > , 


(19b) 


(or vz = v~+iv~ if we prefer to think in terms of two real, 
K k k k 

linearly independent mode-functions '(+ and i+.) Here 

we have defined z = z/v6Mp\. Also, the functions v^(rf) 
have been normalized so that v^v-)' — u^(vi) = i, with 
the convention that in the asymptotic limit H/m —> 0, 
k/(am) —> 00 , we recover standard massless plane waves, 
i.e. vAyj) = e~ lkr] I\[2k in Eq. (I19al) . In the above ex¬ 
pression A(k) is a phase factor, |A(fc )| 2 = 1, and B(k) 
is a real number necessary to connect the two regimes in 
Eqs. (I19al) and (119bl) at k ~ am. For our purposes in this 
paper we will not need to determine these quantities. 

Notice that, while factors of aH/k are small in 


Eq. (HMD, this is not necessarily true for factors of k/aH 
in the case of Eq. (119bl) . where they can even dominate 
for modes inside the Hubble radius, i.e. those modes with 
aH < k < am. This implies that, for modes larger than 
the Compton wavelength, the lowest nonvanishing con¬ 
tribution to the Mukhanov-Sasaki variable is not given 
by limfc^o v k(v) ~ z — iz J d(mr])/z 2 , as one could have 
naively expected. Incidentally the extra terms will be 
crucial to determine the behavior of the modes relevant 
for structure formation. Since aH ~ 77“ 1 and am ~ r] 2 , 
eventually all modes reach this regime. This completes 
our discussion of the Mukhanov-Sasaki variable. 

However, the Mukhanov-Sasaki variable is only an aux- 
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iliary field (remember that this quantity is related to 
the perturbation in the scalar field evaluated in the spa¬ 
tially flat gauge). In order to make contact with obser¬ 
vations we need to move our attention, for instance, to 


the Newtonian potential d', or to the contrast in the en¬ 
ergy density, Se/eq. Introducing the expressions for the 
Mukhanov-Sasaki variable, Eqs. on, into e.g. Eq. (114a[) . 
we obtain for the Newtonian potential 4 


*«<”> = Vl 3S * 1 


i (aH\ 1 


V2k \ k J a 


m 


cos(mf) e~ lkr) , 


-B(k)+iB~ 1 (k) 

5 


faH\ 

/ am \ 

i ■ 

{ k ) 

virJ 

a 3 


if k^> am , 
if k -C am . 


, \j2m 

whereas from the Hamiltonian constraint (112al) . we find for the contrast in the energy density 5 


8e 

£o 


2 1 


3 Mpi 


x 


. y/2m 


—J= ( - cos (mt) e lk7] , 

y/2k \ aH) a 

m 

y/2 m 

A(k) 




-~B(k)+9iB~ 1 (k) 

5 


(aH\ 

/ am \ 

i ■ 

V k ) 

VY J 

a 3 . 


if k am , 
if aH -C k -C am , 
cos 2 (mi) , if k -C aH . 


( 20 ) 


( 21 ) 


Note that the two expressions above have been re¬ 
ported only to the lowest non-vanishing order in the 
series expansions. A pattern of small-amplitude, high- 
frequency oscillations are expected around the expres¬ 
sions in Eqs. m and m, but this will be enough for 
the purposes of this paper. 

For modes larger than the Compton wavelength, 
k < am , the Newtonian potential mimics the behavior 
of standard CDM: one of the solutions remains constant, 
while the other decreases in time with H/a ~ rj~ 5 ; see for 
instance Eq. (7.53) in Ref. [l| for details. On the other 
hand, for modes smaller than the Compton wavelength, 
k > am, the solution oscillates in time with a decreasing 
amplitude, H ~ 1 /rf. This decay in the amplitude of 
the Newtonian potential is characteristic of a barotropic 
perfect fluid p = p(s) -C £ with a nonvanishing Jeans 
scale, c 2 s ^ 0. 

Something similar happens for the contrast in the en¬ 
ergy density, where high-frequency modes with k > am 
oscillate in time with damped amplitude, l/a 2 H ~ t? - 1 . 
In contrast, there appear two different regimes for the 
modes larger than the Compton wavelength. Perturba¬ 
tions larger than the Compton length but still smaller 
than the Hubble radius, aH < k < am, can grow in 
time as 1 /(aH) 2 ~ rf, or decrease as 1 /(Ha 3 ) ~ r]~ 3 , 


Here we have used the identity 


%(*?) = ~ SnG 


n 

2 k 2 


( 22 ) 


that relates, in Fourier space, the Mukhanov-Sasaki variable to 
the Newtonian potential. 


which is again the same behavior as in the standard 
CDM scenario. On the other hand, those modes larger 
than the Hubble radius, k < aH, freeze with a time- 
dependent modulation in cos 2 (mt). It is possible to trace 
back this oscillatory dependency to a background term, 
%' — H 2 ~ cos 2 (mt) / a, that appears in the expression for 
the contrast in the energy density, see Eqs. (E3l) and EH) 
in footnote [5] These oscillations are therefore not ex¬ 
pected to appear in the case of a complex scalar field for 
which H' — H 2 ~ 1/a. 

In particular, the growing mode in the energy density 
contrast 8e/e o ~ 1 /(aH) 2 ~ rj 2 ~ f 2 / 3 (i.e. Se/e o ~ a to 
the lowest order in the expansion series) at those scales 
larger than the Compton wavelength but smaller than 
the Hubble radius, 1/m < A < 1 /H, is usually identified 
with structure formation in the universe (see for instance 
Eq. (7.56) in Ref. Q)- Note however that the (large) os¬ 
cillations for the contrast in the energy density of modes 
larger than the Hubble radius are not present in the case 
of CDM. 

At this point it is interesting to stress once again the 


5 Here we have used the identity 



h ’- h 2 \ i /vgy 

k 2 ) H V z ) 



(23) 


that relates, in Fourier space, the Mukhanov-Sasaki variable to 
the contrast in the energy density. Using the evolution for the 
background universe we can write 


1 

3 


n 1 -h 2 

k 2 


i 

3 




2 

cos 2 (mt) + 0(H/m) 


(24) 
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FIG. 3: Evolution of the contrast in the energy density, 
&/£o|£, as a function of conformal time mr] for three dif¬ 
ferent Fourier modes, k/m = (10 2 , 1/2, 10 -3 ) x IO 273 , in the 
universe depicted in Figure(2] The lower and higher wavenum¬ 
bers were chosen in order to lie outside of the region [aH, am] 
during the evolution. Solid lines represent the results of the 
numerical evolution, whereas dashed lines are obtained from 
an appropriate combination of the envelopes of the linearly 
independent solutions in Eq. m • Modes shorter than the 
Compton wavelength (top panel) oscillate with damped am¬ 
plitude. On the contrary, modes larger than the cosmological 
horizon (bottom panel) freeze, showing an oscillatory behav¬ 
ior associated to the inherent oscillation of the scalar field. 
Finally, those modes that are larger than the Compton wave¬ 
length but smaller than the Hubble radius (middle panel) 
grow at the same rate as the scale factor, as it is expected in 
the case of CDM. These modes can give rise to structures in 
the late universe. Note that since the equations are linear we 
can rescale the vertical axes in the figures arbitrarily, as long 
as the amplitudes remain always less than unity. 


and then in this gauge the contrast in the energy den¬ 
sity grows with time as the scale factor, even for those 
modes that are well outside the Hubble horizon. (It is 
in this gauge that the matter power spectrum is usually 
presented in the literature.) Note that in this coordinate 
system the large oscillations in the energy density dis¬ 
appear, and then it is not clear for us if they could be 
observable in practice or if they are only a gauge artifact. 

Figure [3] shows the evolution of the contrast in the en¬ 
ergy density, 5e/e o|g, for three individual Fourier modes 
that are representative of the different regimes. For this 
figure we have considered the same background evolu¬ 
tion of Figure 0 so that initially we have mr) = 10 1//3 , 
a = 10 2 / 3 and 'H/m = 2 x lO -1 / 3 . In order to calculate 
the contrast in the energy density we integrate numeri¬ 
cally Eqs. (114bl) and ijTHll for the function v £ and the New¬ 
tonian potential H/g, simultaneously with the background 
evolution. We then use the Hamiltonian constraint (1 12 all 
to find the perturbation in the energy density, 5er. 

The top panel in Figure [3] corresponds to a mode 
with k/m = 10 8 / 3 ss 464, which is clearly in the regime 
k > am for the time interval mrj £ (2,11). We can see 
that the contrast in the energy density oscillates with an 
amplitude that decays like l/(a 2 H), as expected. One 
can also clearly see from the figure that we have oscil¬ 
lations with two quite different frequencies: a high fre¬ 
quency oscillation coming from the term e~ lkri , modu¬ 
lated by a lower frequency oscillation that corresponds 
to the term cos (mt) in Eq. (1311) . The middle panel corre¬ 
sponds to a mode with k/m = 1/2 x 10 2 / 3 ~ 2.32, so that 
aH < k < am , also in the same time interval. In this 
case we see that the contrast in the energy density grows 
as 1 /(a 2 H 2 ), with small oscillations. These modes are 
interesting for structure formation since they grow as if 
they were made of CDM. Finally, the bottom panel cor¬ 
responds to a mode with k/m = lO -7 / 3 ~ 4.64 x 10 —3 , 
for which we have k < aH in the evolution. In this 
case the contrast in the energy density rapidly reaches a 
regime where it oscillates with constant amplitude. No¬ 
tice that in all these three scenarios the frequency of the 
oscillations increases as time goes by. This is again to be 
expected since the frequency should be constant in cos¬ 
mological time t, so that it increases in conformal time ry. 


fact that these statements are only valid in the conformal- 
Newtonian coordinate system, where the gauge-invariant 
fields Ik and Se/e o take the same values as the Newto¬ 
nian potential and contrast in the energy density, respec¬ 
tively. If we move to e.g. the synchronous coordinate 
system <j> = B = 0, the contrast in the energy den¬ 
sity is related to the gauge-invariant field de/eo through 
5e s /e o = Se/eo — e' 0 /{a£o) J a’Fchy, the constant of inte¬ 
gration in this formula corresponding to an unphysical 
fictitious mode. As in the standard CDM scenario, the 
expression for intermediate wavelengths aH <C k <C am 
in Eq. m describes now all the modes larger than the 
Compton wavelength of the scalar particle, k -C am , 


IV. DISCUSSION 

For a perfect fluid with equation of state p = p(e) <C e, 
there is a close relation between the Hubble radius and 
the Jeans length, A j ~ c s H~ 1 . This is no longer true in 
the case of a canonical scalar held, where the character¬ 
istic Jeans length is fixed instead by A j ~ c s a\z" / z\~ x ! 2 . 
When this length scale is much smaller than the Hubble 
radius, perturbations can grow. For a perfect fluid this 
is possible only if c 2 -C 1. For a scalar held we have 
c 2 = 1, but we could still satisfy a\z" / z] -1 ^ 2 <C iT _1 . 
This is what happens, for instance, when the zero-mode 
of the scalar held is oscillating with high frequency (when 

















compared to the expansion rate of the universe) around 
a minimum of the potential, where a\z" / z^ 1 ^ 2 3 4 5 ~ m -1 . 
Note that this value for the Jeans length at the scale of 
the Compton wavelength of the scalar particle cannot be 
resolved in terms of a nonrelativistic analysis, as it was 
previously done by Hu et al. in Ref. 5j. Indeed, the 
value of the Jeans length we identify in this paper does 
not coincide with the naive estimation they reported in 
Eq. (4) of that reference, and which has been frequently 
used after that (see e.g. Ref. [To|). 

For those modes larger than the Jeans length the scalar 
field follows the evolution in the standard CDM sce¬ 
nario, except for the large oscillations of the contrast in 
the energy density for modes larger than the Hubble ra¬ 
dius when evaluated in e.g. the conformal-Newtonian 
gauge. Note however that these (large) oscillations are 
not present in the behavior of the Newtonian potential, 
for which the scalar field behaves like CDM, or even in 
the super-Hubble modes of the contrast to the energy 
density when evaluated in e.g. the synchronous gauge. 
As far as we know these oscillations of the contrast in 
the energy density for long wavelength modes in certain 
coordinate systems had not been previously reported in 
the literature, and it will be very interesting to explore 
if they are only a gauge artifact or if on the contrary 
they could have imprinted some signatures on the cos¬ 
mological observables. We expect these oscillations will 
not affect the large scale structure of the universe. Af¬ 
ter all, the distribution of galaxies is only sensitive to 
the subhorizon modes, and these modes follow the stan¬ 
dard CDM evolution. However, they could affect e.g. the 
cosmic microwave background photons at large angular 
scales, which traced super-Hubble scales for a long pe¬ 
riod of time during the cosmological evolution. We leave 
a more detailed analysis of this point for a future work. 

For those modes smaller than the Jeans length the evo¬ 
lution cannot bring the small perturbations in the early 
universe to the nonlinear regime, and the inhomogencities 
are erased. This will introduce a cutoff in the mass power 
spectrum for the distribution of galaxies in the universe. 
Something similar happens in warm DM scenarios fl5j . 


If the mass of the scalar particle lies at the scale of 10 /reV 
or above, as it is expected for the QCD axion [l(|, the 
Jeans length would be smaller than a centimeter, and 
the growth of cosmic structures would be probably in¬ 
distinguishable to that in the standard CDM scenario, 
at least while in the linear regime (see Ref. [§] for the 
case when nonlinearities become important). However, 
if we consider ultralight scalar particles of masses as low 
as 10 -22 eV dUzl, the Jeans length grows to the scale 
of parsecs. This could have observable physical conse¬ 
quences in cosmology [1(1 EH) alleviating, for instance, 
the missing satellite discrepancy [l9| . 

In order to determine properly the new expression for 
the mass power spectrum we would need a more elaborate 
analysis that includes, on the one hand, a knowledge of 
the initial conditions for the scalar field after inflation, 
as well as the evolution of the perturbations during the 
radiation dominated era. We leave this study for a future 
paper. 

Note added .—After the first version of this paper was 
submitted for publication we became aware of Ref. J20| , 
where the linearized Einstein-Klcin-Gordon system is 
considered in the context of cosmological reheating. Al¬ 
though there are some similarities in the analysis, the mo¬ 
tivation of our paper is different. We thanks Prof. James 
P. Zibin for pointing this reference out. 
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